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Abstract 

The annihilation of massive scalar particles in one photon on de Sitter expanding universe is 
studied, using perturbation theory. The amplitude and probability corresponding to this process 
is computed using the exact solutions of the Klein-Gordon and Maxwell equations on de Sitter 
geometry. Our results show that the expression of the total probability of photon emission is a 
function dependent on the ratio mass/expansion factor. We perform a graphical study of the 
total probability in terms of the parameter mass/expansionfactor, showing that this effect is 
significant only in strong gravitational fields. We also obtain that the total probability for this 
process vanishes in the Minkowski limit. 
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I. INTRODUCTION 


The theory of quantum fields on curved space-times has lead to many important results. 
Starting with the paper of Schrodinger from 1939 related to the production of scalar particles 
as a result of the expansion of the background, many authors have turned their attention to 
this subject, using different approaches (l| [ 2 ] [ 6 ] | 8 ]- 17]. One of the most studied subjects is 
re production of scalar particles on Friedmann-Robertson-Walker (FRW) space- 


relatec. 

times 


o t 


1] [15] 17]. The studies of Parker based on minimally coupled field equations, showed 


that particles are produced in pairs and this process is significant only in the strong gravita¬ 
tional fields of the early Universe 0-0 . A similar result was obtained in Ref. [l^] where two 
different mechanisms of particle creation were compared: the production of electromagnet- 
ically interacting particles and free particle creation i.e nonconformally coupled scalar field 


with gravity. This study led to the conclusion that the first type o: 
based on perturbation theory, is dominant in the early Universe 


particle creation, which is 


17]. We also mention that 


a detailed study of fermion pair production from a photon and the time reversed process 


in the FRW metric with conformal flat line element can be found in Ref. 


10]. In the more 


recent papers which study the production of particles on de Sitter space-time, based on the 
S matrix approach | 8 ] 16], the results also prove that the effect of particle production is 
important only for large expansion factor corresponding to the early Universe. 

This letter is focused in the study of the scalar pair annihilation into a photon <p + ip* —» 
7 on de Sitter geometry. The time reversed process is obtained when the scalar pair is 
generated by a photon. We have denoted by p the scalar particle, tp* the scalar antiparticle 


and 7 the photon. It is known that, the QED formalism from 
adapted for the study of QED processes in curved spacetimes 


Minkowski space can be 
Then the definition 


18 


of the transition amplitude in de Sitter geometry can be established [ 6 ]-| 8 | 


16]. In our study, 


the in and out fields are exact solutions of the Klein-Gordon and Maxwell equations in de 
Sitter space, written in the momentum-hclicity basis. 

In section II we present the main steps for computing the expression of the probability 
density for scalar pair annihilation in one photon. Section III is dedicated to the study of 
the total probability of the created photons. This is one of the few studies in which the total 
probability can be calculated exactly. Further we graphically analyse the behaviour of the 
total probability of photon emission as function of the parameter m/u, and given values for 
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the momenta of the scalar particles. In section IV we present the conclusions. 


II. THE TRANSITION AMPLITUDE 


The expression of the line element which describes the de Sitter universe is 19]: 

1 


ds z = dt z — e 2ut dx z = 


(u)t c 


■( dt c — dx ), 


( 1 ) 


t c E ( oo, 0). 


where tu > 0 is the expansion factor and the conformal time is t c = — 

In order to study the annihilation of the scalar pair into a photon in de Sitter geometry, 
we will use the solution of the Klein-Gordon equation from Ref. [5]. Therefore the positive 
frequency solution of the Klein-Gordon equation, in momentum basis is {f|: 


U 0 ) = 


1 [W e-^V m (1) /p A 
2V w (2 tt)3/= “G ) 


: ip X 


( 2 ) 


where H-! is the Hankel function of the first kind and k = \ x 2 — f, where x — — > |. 

k V 4 n n u 2 n 

The theory of electromagnetic field in de Sitter metric was studied in Refs [3f|-14|. In Ref.|3| is 
constructed the theory of the free Maxwell field on de Sitter space-time in the conformal chart 
{t c ,x} and Coulomb gauge i.e. A 0 = 0, {\f— gA l ). }i = 0. Using the conformal invariance, the 
positive frequency solution of the free electromagnetic field in momentum-helicity basis can 
be easily obtained jsj]: 

, , 1 1 


W q,\\ X ) = 


3 ,—iqtc+iqx i 


z l \(q)e 


—2 cot 


( 3 ) 


(2tt)3/2 ^ 

where £\(q) is the polarization vector, which is orthogonal on the momentum of the photon 
q ■ £\(q) = 0, and satisfies the following relations [3]: 


£\(q)z*\'(q) = < 5 aa ' 

= hj- 


qiqj 


( 4 ) 


Further we are interested in studying the annihilation of the scalar pair into a photon using 
the perturbative methods adapted to de Sitter geometry. It is worth mentioning that this 
process is forbidden in the scalar QED from Minkowski space, because of the simultaneous 
conservation of energy and momentum. The transition amplitude for this process is 16]: 


= ~e d 4 x^g(f^(x) di fp>(x))w% x (x), 


( 5 ) 
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where y/^g = e 3ut and f d g — f(dg ) — g{df ) is the bilateral derivative. After we replace 
the solutions of the free scalar field from equation (J2J) and Maxwell field from equation (J3]) 
in the expression (|5|) , the transition amplitude will have the following form: 

(p' -p)e*\(q) f + °° 




ie e 

16 y/n y/q (27t) 3 

r+oo 

x / d 3 x e l( P + p'-^ s . 


dt 


- e iqtc H ( { } (?-e 

L0 lk \0J 


4 1} ( v -*~* 


( 6 ) 


The spatial integral gives a delta Dirac function of the form: (2n) 3 5 3 (p+p'—q), which assures 
that the momentum conservation law is preserved in this process. The temporal integral 
contains the contribution of the expansion of the space in this process, and is expressed in 
terms of Hankel functions. If we consider the new variable z = —f c , this type of integral 
can be solved using equation (H71) from Appendix, which relate the Hankel functions and 
Bessel J functions. The exponential function e~ iqz can be expressed in terms of the modified 
Bessel function K\ (z), using the relation (1TB]) from Appendix. Then the expression of the 
temporal integral becomes: 

r- 7T [ dzz 3/2 Ki(iqz)[J^ ik (pz)J- ik (p'z) + e 2nk J ik (pz)J ik (p'z) 

V7rsmh (nk) J 0 2 

-e nk J-i k {pz)J ik {p'z) - e nk Ji k (pz) J-i k (p'z)\. (7) 


The result of the integrals which contain Bessel J functions with equal indices depend on the 
derivative of the Legendre functions of the second kind Q v (u), with respect to the variable 
u = P + yv~ g ’ as one can see l rom l^ e equation (TT9T) . Further using equation (l2Tj) . the final 
result of these types of integrals can be expressed with hypergeometric Gauss function 2 Fi- 
For solving the integrals which contain Bessel J functions with opposite indices we use the 
relation (1201) . The result for these types of integrals depends on the Appell function F 4 . The 
final expression for the transition amplitude is: 

7 = 6 3 (p + p' -q) £*x(q) 

x [gfc (p, p',q) + g -kip, p',q) + (p,p',q) + h -fc ip, p',q)\, ( 8 ) 


where the functions g ±k (p,p',q) and h ± k {p,p',q) are defined as follows: 


g ± k (p,p',q) 


iyjq (k 2 + |) e ±nk 

4 (pp') 3 / 2 cosh(7T k) sinh 2 (7rfc) 


ie ±nk 2 Fi 


+ 2 A] 


/ 3,3. 1 + u\ 

{2 ±ik ’2 Tik ' 2 '—) 


3 3 

- ±ik ,~-F ik ; 2; 
2 2 


1 -u\ 

—J 


(9) 
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h ±k{p,p',q) = t Q ]( n (^ 7 ) ^4 1 =F - *0, - *oV (10) 

7tk smn(7TK) / \2 g g / 

The probability of transition is obtained summing after the photon polarizations and squar¬ 
ing the transition amplitude (J8]) . If we use the identity \S 3 (p + p' — q ')\ 2 = V5 3 (p + p' — q'), 
the final result of the probability density is: 

P<p+V-h = -p)el(v)\ 2s3 (p + p' - 9 ) 

A A 

[\g k (p,p',q )\ 2 + \g-k(p,p'i<i )\ 2 + I K(p,p',q )\ 2 + \h-k(p,p',q )\ 2 
+ g t (p, p q)g-k (p, p',q) + g Kp, p q)K(p, p',q) + g* k (p, p g)h_*(p, p q) 

+ g-fc(p,p', g)gfc(p, p' , 9) + g*fc(p,p', ?)hfc(p,p', q) + g*- k (p,p' , q)h-k(p,p ', g) 

+ hfc(p,p / ,g)g fe (p,p / ,g) +hfc(p,p',g)g_ fe (p,p / ,g) + hfc(p,p',g)h_ fe (p,p / ,g) 

+ hlfc(p,p / ,g)gfc(p,p / ,g) + bl fc (p,p',g)g_ fc (p,p',g) + h*_ fc (p,p / ,g)h fc (p,p / ,g)] . 

( 11 ) 


We observe that the transition probability is proportional with a term of the form: 

Sa \(p' ~ P)e* x (q )\ 2 . Using the momentum conservation p' = q — p and equation (j3J), this 
term can be expressed as follows fldl ]: 

a / 2 2 

^2\(p' - p)£* x {q)\ 2 = 4p 2 sin 2 6 qp = P f sin 2 9 pp '. (12) 

A ^ 

The final result of the probability of transition depends on the parameter x = m/u and 
the particle momenta p,p',q- This allows us to study the behaviour of the probability in 
terms of m/u for fixed momenta. It is also worth to mention that, using the perturbative 
approach the probability density is a quantity that depends on the momenta of the particles. 
Consequently, a divergence in the integrals over the final momenta could be solved using 
regularization methods. 


III. THE TOTAL PROBABILITY 


The total probability is obtained integrating the transition probability (ITT]) , after the 
final momenta of the photon: 


Plot = / d 3 q P 


ip+ip*^ 7 - 


(13) 
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In the case of scalar pair annihilation, the expression of the total probability can be evaluated 
analytically, because we have to perform only one integral over the final momentum of the 
photon. Taking into account that the result of the probability of transition is proportional 
with delta-Dirac function S 3 (p+p' — q), the momentum integration can be easily performed. 
Then the final result for the total probability per unit of volume reads: 

Plot 


[|gfc(P,p')| 2 + Ig-fcfep')! 2 + |h*(p,p') I 2 + I h_fe(p,p')P 


e 2 sin 2 0 pp i p 2 p' 2 
V 1287T \p + p'\ 2 

+ gk(PiP')g-k(PiP') +gUP,P')hk(p,p') +gl(p,p')h-k(p,p') 

+ g*-k(p,p')gk(P,p') +g-k(p,p')hk(p,p') +g*_ k (p,p')h-k(p,p') 

+ K(p,p')gk(p,p') + K(ppp') h -k(p,p’) + K(ppp')g-k(p,p') 

+ h*-k(P,P')gk(P,P') + h*_ k (p,p')g_ k (p,p') + h*_ k (p,p')h k (p,p')] . 


(14) 


The expressions of g ± k (p,p') and h ±k(p,p') which are now functions only of p and p' are: 


g ±k(P,P') = 


dz7r/c 


i\p + p'\ 1/2 {k 2 + j)e 

4 (pp') 3 / 2 cosh(7r/c) sinh 2 (7rA:) 


( § ± | =F it; 2;i±^£V 


"3 3 1 

2-^1 ( X ± ik, - =F ik] 2; — 


cos Q 


pp 


(15) 


h ±k(p,p') = 


| p + p 


71-5/2 


7T k sinh(7r/c) \p 


P 




F 4 ( 11,1 ± ik, IT ik; 


iO, 


P 


/ 2 


\P + P'\ 


iO 


(16) 


where \p + p'\ = (p 2 + p’ 2 + 2pp’cosdpp/) 1 ^ 2 . 

Further we want to graphically analyse the total probability of photon emission as a 
function of x — m/cu G (1.5, oo), for different values of the angle 9 pp >, keeping p and p ' 
fixed. As can be seen from (ll5il the hypergeometric function 2 F 1 is divergent for 6 pp i = 
0, 7 r. In this situation we will study the behaviour of the total probability for small angles 
(7 t/25 < 9p P / < 7r/15) as well as for large angles (0.97T > 9 pp i > 0.727r). The h±k functions 
are proportional with a factor {k sinh(Trfc))^ 1 , which make them convergent for large m/uo. 
The relevant contribution at large m/cu will be given by this factor and for that reason we 
approximate the Appell functions F 4 ~ 1. 
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FIG. 1: The total probability P as function of x = m/uj , for p = 1 = 2 and 9 pp ' G {^, §}. 
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FIG. 2: The total probability P as a function of x = m/u, for p = l.p' = 2 and small values of 
the angle 9 pp i. 








-0=O.72n. 0=0.8371 - 0=0.971 


FIG. 3: The total probability P as a function of x = m/ui , for p = l,p' = 2 and large values of 
the angle 9 PP '. 

From Fig.l it can be seen that, in the particular cases when 9 pp t e {f, f, the total 
probability is a function which decreases as the parameter m/uo increases. This property of 
the total probability is also preserved for the small/large angle cases. The obvious conclusion 
is that the total probability of scalar pair annihilation is significative only for m/uj 6 (1.5, 2.5) 
where the gravitational fields are still strong. From Fig.2 we see that, in the case when 9 pp i 
is close to zero, the total probability is large compared to the other studied cases. This 
important observation must be interpreted by analysing the analytical formula (fTbl) of the 
total probability. As it can be seen from (fT6l) the total probability of scalar pair annihilation 
is proportional with sin 2 9 pp i , which arises from the polarization term. One is tempted to say 
that when 9 pp i —> 0, the total probability also approaches zero, but the functions g±k(p,P') 
have a divergent behaviour for 9 pp t = 0, because the Gauss hypergeometric function 2 F\ 
becomes also divergent in this case. Also in the large angle case, the total probability 
of massive scalar pair annihilation increases when 9 pp i —> n, but this probability is small 
comparatively with the small angle case. This result shows that, it is more probable for the 
scalar pair to annihilate into a photon when the momenta of the scalar particles are parallel 
and orientated in the same direction. Also from these graphs we can draw the conclusion 
that the annihilation probability is minimal when 9 PP > = f. Another important observation 
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which can be made from Figs. (1) —(3) is that, the total probability of scalar pair annihilation 
vanishes in the Minkowski limit, when the parameter m/uj —> oo. The Minkowski limit can 
be also verified by analytical calculations. From the expressions (1131) — (fltrl) it can be seen 
that when m/co 1, the functions g k (PiP') and h k (p,p') decreases rapidly to zero with a 
term of the form . 

It seems from our study that the total probability of the created photons is divergent for 
9 PP i = 0,7r. A possible way to eliminate the divergence is by using regularization methods, 
since the density of probability depends on the final momenta. 

IV. CONCLUSION 

In this letter we have studied the annihilation of massive scalar particles into a photon 
in de Sitter spacetime, using a direct perturbative calculation. We were able to obtain an 
analytical formula for the total probability of scalar pair annihilation. The total probability 
of created photons is a function which depends on the parameter m/uj and the momenta of 
the scalar particles. Also, the graphical results show that the total probability is large when 
the scalar particles annihilate at small angles. In the limit uj —)■ 0 we recover the flat space 
result where the probability for this process vanish. 


V. APPENDIX 


For solving our temporal integral we express the Hankel functions in terms of Bessel J 
functions as follows 


20J: 


= 

A 2, W = 


j- p {z) - e infl J p (z) 
i sin Tip, 

e^J^z) - J-p(z) 
i sin Tip 


(17) 


To replace the exponential from our integrals, we use the expression of Bessel K function of 
index 1 |20|: 


aw= v£ e -- 


(18) 


The integrals which contain products of Bessel J functions with equal indices can be ex¬ 
pressed in terms of the Legendre functions of second kind Q u (u), where 2 abu = a 2 + b 2 + c 2 
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20 |: 


dxx 2 Ki(cx) J u (ax) J u (bx) = 


C 2 


d 


Qu~b( u ) 


(19) 


a/ 27 r ( ab ) 3 4 5 6 7 8 9 / 2 du 

The integrals which contain products of two Bessel J functions with opposite indices can be 
solved using the relation 20]: 


dxx^Ki (cx) J v (ax) J_ v (bx) = c 2 ^ 

V 27 tu 


S n 2 h 1 2 

( 20 ) 


where F 4 is the Appcll hypergeometric function with double argument. The above integrals 
are convergent for 7 Ze(c) > 0 and for solving the integrals we take c —> e — iq with e > 0 and 
finally we consider the limit e —> 0. 

The derivative of the Legendre function Q v (u) in respect to u can be written in terms of 
the Gauss hypergeometric function 2 Pi as follows js ]: 


d n , , . n s vrz/(z/ + 1) 

Qv(u i zO) . 
du 4 sm itu 


2 F] l-i/, 2 + 1/; 2; 


u 


+ 2^1 I 1 — u, 2 + 1 /; 2; 


1 + M 


This relation help us to establish the final form of the transition amplitude. 
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